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Abstract 

The electromagnetic properties of the deuteron are investigated within a Light-Front Hamilto- 
nian Dynamics framework, with a current operator that contains both one-body and two-body 
contributions. In this work, we are considering new two-body contributions, with a dynamical na- 
ture generated within a Yukawa model and a structure inspired by a recent analysis of the current 
operator, that acts on the three-dimensional valence component and fulfills the Ward-Takahashi 
identity. Preliminary results for the magnetic moment are shown. 



1 



I. INTRODUCTION 



A careful description of the deuteron, retaining only a finite number of baryonic and 
mesonic degrees of freedom, could help to single out new, experimental signatures, related 
to the underlying degrees of freedom (see [l, 2] for recent reviews of the QCD effects in the 
deuteron). In the present context, careful description means an approach that is able to 
fulfill general properties, like the extended Poincare covariance (the discrete symmetries are 
included), Hermiticity and current conservation j^|. 

Deuteron electromagnetic (em) observables offer a valuable play ground for testing theo- 
retical ideas (see, e.g., |4(). In particular, the challenges for theorists appear very stimulating 
in view of the Th. Jefferson Lab. (TJLAB) upgrading to 12 GeV, that should open very 
intriguing scenarios, as described, e.g., in PAC-34 and PAC-35 Reports [5]. But, the list of 
issues to be coped with is long, if one would like to address non standard effects, e.g. like 
the 6-quark bag Q]. In our opinion (see, also ^|), efforts should be invested on the 
analysis of: i) the consistency between dynamics and operatorial structure of the current; ii) 
the strong interplay between different ingredients, e.g. between the dynamical content of the 
two-body currents and the nucleon form factors; iii) the two-photon exchange effects, that 
could affect the extraction of the em form factors (though recent estimates 6] assign them 
a minor role, ~ 1%); iv) further baryonic degrees of freedom, like isobar configurations, till 
now investigated only within a non relativistic approach [7]. 

Our aim is to include new dynamical two-body contributions into the em current, in a 
Light-Front Hamiltonian Dynamics (LFHD) framework (see, e.g., js| for a review of the Rel- 
ativistic Hamiltonian Dynamics), for describing the deuteron em observables, still satisfying 
;he extended Poincare covariance. This work expands the investigation carried out in Refs. 
9), [loj], where a current operator containing the one-body term and a two-body contribu- 
tion, needed for satisfying the Hermiticity, was considered. In the present approach, we add 
two-body terms with a dynamical nature, generated by the presence of an explicit one-pion 
exchange (see also 111]). These two-body currents are inspired by an exact analysis of the 
four-dimensional (4D) current corresponding to a Yukawa model in ladder approximation, 
where two fermions interact with a pseudoscalar, massive boson 
an approximation, there is no photon-boson coupling. The 4D current corresponding to the 
field theoretical model is projected onto the three-dimensional (3D) LF hyperplane, so that 



121 ] (notice that, within such 
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one obtains an operator, that i) acts on the 3D LF valence component of the interacting- 
system state and ii) automatically fulfills the Ward-Takahashi Identity (WTI). Moreover, 
in Ref. 12], it has been shown that one can properly truncate, in the Fock space, the LF 
current and still be able to satisfy the correspondingly truncated WTI. In particular, the 
Fock expansion is ordered in powers of the interaction. In order to improve the calculations 
of the em deuteron observables of Refs. js), [lo|], in this work we consider the first-order 



LF current operator, obtained by applying the approach of Ref. [12j to an interacting La- 
grangian £ = —igps^^^r^ ■ 0, where \l/ is the fermion field and an isovector pseudoscalar 
boson field. 

In Sec. [Ill the choice of the reference frame and some generalities of our approach will 
be presented. In Sec. IHH the current adopted in our LFHD approach will be illustrated. In 
Sec. IIVt the preliminary results, with only a part of the two-body dynamical contributions, 
will be shown for the deuteron magnetic moment. In Sec. IVl summary and perspectives will 
be presented. 

II. CHOOSING THE FRAME 

Our theoretical frame is the LF Hamiltonian Dynamics combined with the Bakamjian- 
Thomas construction [3] of the Poincare generators for an interacting system. Within such 
a framework, as it is well known, a relativistic square mass operator of a two-body system 



can be immediately identified with the Schrodinger equation for that system 14|. This 



allows one to formally embed the standard non relativistic deuteron wave function into a 



fully Poincare covariant description. The on 
part, where the Melosh rotation operators 



y 



particular care is the treatment of the spin 



15] are requested (see, e.g., [10| for details). 



Therefore, within the previous approach, one can rigorously fulfill the Poincare covariance, 
for an interacting system with a finite number of degrees of freedom. Heuristically, one could 
say that the Relativistic Hamiltonian Dynamics with a Bakamjian-Thomas construction 
allow one to implement a description of an interacting system that falls between the non- 
relativistic quantum mechanics and local relativistic field theories. This simple view is 
further strengthened onto the LF hyperplane, since in this case a sharp and clean separation 
between the center of mass motion and the intrinsic one can be straightforwardly achieved. 
Let us also notice that the LF boosts are kinematical transformations and that the spectral 
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condition P + = P° + P z > favors a description with a finite number of degrees of freedom. 

The other ingredient for developing our description of the em deuteron form factors is the 
reference frame where we calculate the theoretical observables. If one had the complete and 
calculable theory at disposal, the choice of the frame should not represent an issue, given 
the full covariance,. Actually one is dealing with an approximate scheme, and therefore 
the choice of the reference frame, where one can more easily implement the constrains for 
preserving the general properties, becomes strategic. Following Ref. {3], one can show that 
in a Breit frame where the momentum transfer is longitudinal, i.e. q± = 0, the symmetry 
of the physical process can be exploited for reducing the constraints imposed on the current 
operator by the extended Poincare covariance to a simple rotational covariance around the 
z = q -axis. Then, in this frame, any operatorial dependence, fulfilling the rotational 
covariance around the z — axis, is allowed in the construction of a current operator, that, 
in particular, i) depends parametrically upon the CM momenta and ii) acts on the intrinsic 
variables. For instance, the matrix elements of the current operator can be approximated 
by the ones obtained from one-body Dirac and Pauli currents, and then, since the rotational 
covariance is safe, it turns out that the extended Poincare covariance can be fulfilled by 
using the proper transformations (Lorentz boost, rotations, translations etc.). Moreover, 
in order to implement the Hermiticity, one has to add a term that contains the dynamical 
generator of the transverse rotations, obtaining necessarily a two-body contribution to the 
current. It should be pointed out that such a contribution must be distinct from the ones 
that we are going to discuss in what follows, since the last ones contain dynamical effects in 
a more explicit way. 



III. PROJECTING THE EM CURRENT ONTO THE LIGHT-FRONT HYPER- 
PLANE 



In this Section, a brief illustration of the results of Ref. 
relevant for the present calculations, is given. 



121 ] . in particular the ones directly 
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In Refs 
LF hyperpiajic (1.0. .r 



17H19|. by using the projection of the 4D physical quantities onto the 3D 
0), and the Quasi- Potential approach for the Transition 



x° + x 3 



Matrix 16], it has been established a formally exact, one-to-one correspondence between i) 
the 4D Bethe-Salpeter amplitude of an interacting system and the 3D LF valence component, 
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ii) the matrix elements of the 4D em current and of the LF current, both fulfilling the proper 
Ward-Takahashi identity. The previous two ingredients allow one to keep separate the trivial 
propagation in the global time (J2i,n x f f° r a system with N constituents), and to focus on 
the dependence upon the relative-time propagation, governed by the internal dynamics of 
the system. In particular, within this approach, the valence component is found to be the 
eigensolution of a 3D dynamical equation, with an effective 3D interaction that is exactly 
related to the 4D kernel of the BS equation. The exact correspondence between 3D and 4D 
quantities is accomplished through operators that contain the 4D interaction and cannot be 
determined in a simple way. This difficulty can be overcome in a workable way by developing 
a possible approximation scheme, for constructing solutions, based on the Fock expansion 
of the relevant quantities. If one truncates the Fock space, one can use the corresponding 
truncated Fock basis, in order to expand quantities, like the effective interaction and the 
current operator, onto the LF hyperplane. Clearly, by using a truncated space only the 
kinematical symmetries can be satisfied within a field theoretical approach (given the infinite 
number of degrees of freedom), but if one restricts to a Relativistic Hamiltonian Dynamics 
approach, where a finite number of degrees of freedom is taken into account, then one can 
adopt the truncated operators and recover the full Poincare covariance. In particular, if 
the truncated current contains operators that satisfy the covariance around the z-axis, and 
furthermore one uses the valence wave functions, eigensolution of the properly truncated 
mass equation, in the evaluation of the matrix elements, then both the current conservation 
and the Poincare covariance can be implemented. It should be stressed that the truncated 
mass operator has to properly commute with the Poincare generators, as requested by the 
Bakamjian-Thomas construction. 

In the actual calculations, there is another relevant issue for the LF projection, that makes 
the application of the procedure to the fermionic case sharply different from the bosonic one. 
The Dirac propagator can be separated in an on-shell term and an instantaneous (in LF 
time !) propagation, viz 

i c (k ) = ft + m = hn + m 7 + m 

K> k 2 -m 2 + te k+(k~ - k~ n + 2k+ V) 

where k^ = (|k^| 2 + m 2 )/k + is the minus-component of k^ n) such that k on ■ k on = m 2 and 
the second term is the instantaneous one, as shown by the Fourier transform in k~ and kj_, 
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i.e. 



J dk~ dk ± exp[-i(k-x + /2 - ik ± ■ xjj] = (2tt) 3 5(x + )5(x±) . 

The instantaneous term has a great impact on the analysis of the matrix elements of 
the current operator, and it produces very peculiar operatorial structures in the many-body 
contribution to the LF current. It should be emphasized that this fact is related to our 
choice of the reference frame, where q + 7^ 0, since in this case the extraction of the em 
form factors involves matrix elements of both the plus and perp components of the current. 
In the frame with q + ^ the Poincare covariance can be demonstrated within a LFHD 



approach [3j, and therefore the famous angular condition issue (see, e.g. |20|) does not any 
more plague the matrix elements. 

In Ref. 12J, the explicit expressions for many-body terms of the em LF current have 
been obtained, within a Yukawa model in ladder approximation for two fermions. In our 
approach for evaluating the em deuteron observables, we take into account the current only 
up to the first-order in the interaction (see the next Section), as shown by the diagrammatic 
representation of Fig. 1. It is worth noting that one has at least three particles in flight (the 
LF-time flows from right to left). 

IV. PRELIMINARY RESULTS FOR THE DEUTERON MAGNETIC MOMENT 

In a Breit frame where q± = and therefore q + ^ 0, the matrix elements of the em 
current operator for an interacting system can be defined in terms of the LF current and the 
dynamical transverse component of the Bakamjian-Thomas rotation generator S± as follows 

f(qe z ) = + L^M-tt)] e^s.W&T^s. (2) 

where L£[r x ( — 7r)] is the 4D representation of the 7r-rotation around the x-axis, and generates 
a term necessary for implementing Hermiticity. In Eq. (T5]), the operator J^{qe z ) is given by 

J^z) = 3Le + jlo = n Jo^o + n \VAojf + XA v] n (3) 



where LTo is the so called free reverse LF projector of Ref. [121 ] . that singles out the positive- 
energy sector (modulo some kinematical factors) and Ao = Gf ree — G g i Q b, with Gf ree the 
standard two-fermion free propagator and G g i Q b an auxiliary Green's function, that takes 
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FIG. 1: Diagrammatic analysis of the LF first-order em current, obtained within a Yukawa model 
in ladder approximation. In Ref. [12( one can find more details and explicit expressions. 

into account t he g lobal-time propagation. The operator V is the interaction, mediated by a 



pion (see also 



11 



]) 



V = i g 2 fx ■ r 2 



7l g 72 F [(p2-piY 
[(P2 - Pi) 2 -ml + ie] 



(4) 



with J r [ 
of Ref. 



]k— Pi) 2 } the 4D vertex form factor counterpart of the nonrelativistic vertex function 



2l|. The quantities with hats represent proper operators. In Eq. ([3]), the one-body 



contribution (see ]9 



lOj for details) is obtained from 



i=l,2 L 



(5) 
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where the free nucleoli current is 

r N = -F 2N \{P' -P) 



J 2M 



1JV 



p) 2 ] + i^v[(p'-p) 2 ]) 



(6) 



with Fitv and .F2.1v the Dirac and Pauli form factors, respectively. 

What about current conservation and charge normalization? In the chosen Breit frame, 
current conservation and charge normalization read respectively as follows 



and 



[P f ,d\j+(qe z )\d;Pj = (P f ,d\j-(qe t )\d; P) 



(Pi,d\j + (0))\d',Pd = (Pi,d\± [j + (0) + J-(0)} \d-Pi) = e 



(7) 



(8) 



If WTI is fulfilled, then one obtains the current conservation, once matrix elements are taken 
between eigensolutions of the mass equation constructed from the proper Green' function. 
This is not the case in our phenomenological calculations, since we are adopting the deuteron 
wave functions corresponding to realistic interaction, like CD-Bonn [2lJ| or AV18 22]. But, 
in the elastic processes like the one we are considering, current conservation follows after 
implementing Hermiticity |3|, given by (note the change of the z-axis) 



(P f ,d\f{qe z )\d;Pj = {P h d\f{-qe z )\d-P f Y 



Notice that the charge normalization can be fulfilled if in Eq. (T5]) one defines 



(Pi,d\J-(0)\d;Pj = (P^dU+mdlP] 



(9) 



(10) 



This leads to assume that such an equality holds for any momentum transfer. It should be 
pointed out that for evaluating the em observables only j + (qe z ) and j 1 ^(qe z ) are relevant. 

Preliminary results for the magnetic moment of the deuteron have been obtained by 
retaining only the two-body interaction terms, corresponding to the diagrams depicted in 



the first line of Fig. 1, and with the explicit expressions given in Ref. 



12] (modulo the 



isospin dependence and the pseudoscalar coupling). The nucleon form factors adopted in 



the calculations are the ones of Ref. 



23] . that are in nice agreement with the most recent 



measurements of the proton form factors (see, e.g., 25]) and represent a first microscopical 
interpretation of the possible zero in the ratio n P G p E /G p M for Q 2 > 8 (GeV/c) 2 in terms 
of interference between the valence and non valence component of the proton state. It is 
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worth noting that the behavior of the nucleon form factors, in a wide kinematical region, 
enters in the evaluation of the static (Q 2 — >■ 0) em properties of the deuteron, according 

nn 

to the approach of Refs. [9J, [HJ. In Table 1, the preliminary values, for different deuteron 
wave functions, correspon ding to three realistic NN interactions, are shown. In particular, 
the CD-Bonn |2l| , RSC93 



24J and AV18 22j interactions have been used. 



'he calculations 



corresponding to the one-body contribution with the Hermiticity term [9_|, [10] have been also 
shown. It should be mentioned that, within the last approximation scheme, the experimental 
quadrupole moment is fairly well described, with an underestimate of the order of 4%. 



Interaction 


Pd 




LFD 

r'one 


til? 


CD-Bonn 


4.83 


0.8523 


0.8670 


0.863 ±0.002 


RSC93 


5.70 


0.8473 


0.8637 


0.861 ±0.002 


Avl8 


5.76 


0.8470 


0.8635 


0.860 ±0.002 



Exp. 0.857406(1) 



The results appear encouraging. But, for a vanishing momentum transfer (see the analogous 
discussion for the one-body case in Refs. 0,0), it is necessary an accurate study of the 
numerical convergence of the multifold integrals, that enter the calculations. From the charge 
normalization, one can obtain the probabilities of the valence and non valence components. 
At the present stage we have obtained Prob NV ~ 0.01. 



V. SUMMARY & PERSPECTIVES 



In this contribution, we have presented preliminary results for the magnetic moment of 
the deuteron, including two-body, dynamical contributions [l^] to em current, within the 
LFHD approach proposed in Refs. [3|, |9j, llQj • The approach is fully Poincare covariant. 

The new two-body terms have been inspired by an exact analysis of a Yukawa model, 
in ladder approximation, for two interacting fermions, carried out in Ref. 12j . It turns 
out that, onto the LF hyperplane, one obtains a LF current fulfilling the Ward-Takahashi 
identity, at each order in the Fock expansion. 

The systematic analysis of the deuteron em form factors is started, and it represents a 
non trivial task from the numerical point of view, given the many, multidimensional integrals 
to be performed with very high accuracy, in particular at low Q 2 . First results obtained 
by using i) the deuteron wave functions corresponding to CD-Bonn, RSC93 and AV18 NN 



potentials and ii) the interaction contribution to the two-body current, as depicted in the 
first line of Fig. 1, appear consistent with the expectations, namely a very low probability 
for the component beyond the valence one, and a magnetic moment in fair agreement with 
the experimental values. Calculations of the quadrupole moments, more delicate from the 
numerical point of view, are in progress. 
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